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Abstract 

The leptoproduction of J/ip mesons is studied in inelastic reactions for four momentum 
transfers 2 < Q 2 < 100 GeV 2 . The data were taken with the HI detector at the electron 
proton collider HERA and correspond to an integrated luminosity of 77 pb -1 . Single dif- 
ferential and double differential cross sections are measured with increased precision com- 
pared with previous analyses. New leading order calculations within the non-relativistic 
QCD factorisation approach including colour octet and colour singlet contributions are 
compared with the data and are found to give a reasonable description of most distribu- 
tions. An exception is the shape of the distribution in the J/ip fractional energy, z, which 
deviates significantly from that of the data. Comparisons with photoproduction are made 
and the polarisation of the produced J/ip meson is analysed. 
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1 Introduction 



Inelastic leptoproduction of J/ijj mesons at HERA, e+p — > e+J/ip+X, is, dominated by boson 
gluon fusion, 7*g — > cc. The aim of current experimental and theoretical efforts is a detailed 
understanding of this production process. Before HERA started operation, the limited amount 
of lepto- and photoproduction data ( [[[]] and references therein) was found to be described by 
the Colour Singlet Model (CSM) [g]. In the CSM the cc pair is produced in the hard 7*3 
interaction in the quantum state of the J/ip meson, i.e. in a colour singlet state with spin 1 
and no orbital angular momentum. This is possible due to the emission of an additional hard 
gluon (see Fig. [j]b). The process was advocated as a means of determining the gluon density in 
the proton, since it is calculable in perturbative Quantum Chromodynamics (pQCD) using e.g. 
potential models for the formation of the J / ip meson. 

In recent years the interest in inelastic J / ip production has shifted considerably since the CSM 
fails to reproduce the production rate of J/ip and ip(2S) mesons in pp collisions by more than 
an order of magnitude [|J]. Nowadays, one of the main aims is the investigation of the role of 
"colour octet" contributions, which have been invoked to describe the pp data. Colour octet 
contributions arise naturally in the theoretical description of quarkonium production based on 
non-relativistic QCD and factorisation (NRQCD) [Q]. NRQCD is an effective field theory in 
which the J/ip production process factorises into terms for the short distance transition (e.g. 
7*g — > cc(g)) and long distance matrix elements (LDMEs) for the transition of the cc pair into 
an observable meson. The cc pairs can now be in many different angular momentum states, in 
colour singlet and also in colour octet states, in which case the transition to the J/ip meson is 
thought to proceed via soft gluon emission. The short distance coefficients are calculable in 
pQCD and a double expansion in the strong coupling parameter a s and v, the relative velocity 
of the quark and antiquark, is obtained. Many contributions are possible (examples are shown 
in Fig. P and only the most important contributions are kept in a specific calculation 
The leading term in the velocity expansion is the colour singlet term, so if it is assumed that 
all other terms do not contribute, the CSM is recovered. Although the octet LDMEs are at 
present not calculable, they are assumed to be universal. They have been extracted from the 
measurement of J/ip production in pp collisions by fitting the leading order (LO) theoretical 
calculation to the data (see e.g. [0, ||| and references therein) and are then used in predictions 
for electroproduction. 

First attempts to establish the relative importance of colour octet contributions in lepton proton 
interactions were made in the photoproduction limit, Q 2 — > P,[T0t], where Q 2 is the negative 



squared four momentum transfer. The predicted large contributions at high values of the J / ip 
fractional energy, z, were not observed. Here, z denotes the J/ip energy relative to the photon 
energy in the proton rest system. In the analysis of data at high Q 2 the dependence of the cross 
section on Q 2 may give additional insight into the production process [j^]. 

Analysing leptoproduction at finite Q 2 has experimental and theoretical advantages compared 
with photoproduction. At high Q 2 theoretical uncertainties in the models decrease and resolved 
photon processes are expected to be negligible. Furthermore background from diffractive pro- 
duction of charmonia is expected to decrease faster with Q 2 than the inelastic process. The dis- 
tinct signature of the scattered lepton makes the process easier to detect. A first comparison be- 
tween data and NRQCD calculations in the range 2 < Q 2 < 80 GeV 2 and 40 < W < 180 GeV 
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Figure 1: Generic diagrams for charmonium production mechanisms: a) Photon gluon fusion 
via a "2 — > 1" process; b-d) "2 — > 2" processes, a-d) contribute via colour octet mechanisms, 
while b) can also contribute in colour singlet processes. Additional soft gluons emitted during 
the hadronisation process are not shown. 



was presented in [|TT|], W being the mass of the hadronic final state or equivalently the centre 
of mass energy of the photon proton system. The NRQCD calculations compared with the data 
in JTT| ] were performed taking into account only "2 — > 1" diagrams [|5|] (compare Fig. and 
disagreement between data and theory was observed both in the absolute values of the cross 
sections and in their shapes as functions of the variables studied. 

In this paper, an analysis of e + p — > e + J/ip + X is presented in the kinematic region 
2 < Q 2 < 100 GeV 2 and 50 < W < 225 GeV with increased statistics compared to our 
previous publication [|TT|]. Differential cross sections are measured for the whole Q 2 range and 
for a subset with Q 2 > 12 GeV 2 . The data are compared with theoretical predictions [j^] in the 
NRQCD framework taking into account colour octet (CO) and colour singlet (CS) contributions. 
In contrast to the previous NRQCD calculation, diagrams of the type "2 — > 2" are taken into 
account (e.g. diagrams |l]b, c and d). The J/ip polarisation is measured by analysing the decay 
angular distribution and its Q 2 dependence is investigated. The polarisation measurements are 



compared with the prediction of a calculation [ |12| ] within a "k t factorisation" approach, i.e. 
allowing transverse momentum ("k") for the incoming gluon, using unintegrated parton density 
functions and off- shell matrix elements including colour octet and colour singlet contributions. 
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2 Detector, Kinematics and Simulations 



The data presented were collected in the years 1997-2000 and correspond to a total integrated 
luminosity of 77.0 ± 1.2 ptT 1 . HERA was operated for most of this time with 27.5 GeV 
positrons. Roughly 12% of the data were taken with electrons of the same energy. In 1997 
the proton energy was 820 GeV. It was increased to 920 GeV thereafter (sample of ~ 63 pb _1 ). 

The experimental methods are similar to those described in [JTT]] and further details can be found 
in J/ip mesons are detected via the decays J/ip — > and J/ip — ► e + e~ (branching 

fractions of 5.88 ± 0.10% and 5.93 ± 0.10%, respectively 0). 



2.1 Detector 



A detailed description of the HI detector can be found elsewhere [ p3| ] . Here we give an overview 
of the most important components for the present analysis. The central tracking detector (CTD) 
of HI consists mainly of two coaxial cylindrical drift chambers for the measurement of charged 
particles and their momenta transverse to the beam direction and two polygonal drift chambers 
for measurement of the z coordinates^]. The CTD is situated inside the solenoidal magnet, which 
generates a field of 1.15 T. The tracking system is complemented in the forward direction by 
a set of drift chambers with wires perpendicular to the beam direction which allow particle 
detection for polar angles 9^7°. Multiwire proportional chambers are used for triggering 
purposes. 

In the Q 2 range studied here, the scattered lepton is identified through its energy deposition 
in the backward electromagnetic calorimeter SpaCal [|T6|]. The SpaCal signal is also used to 
trigger the events, in conjunction with signals from the tracking chambers. A drift chamber 
(BDC) in front of the SpaCal is used in combination with the interaction vertex to reconstruct 
the polar angle 9 e of the scattered lepton. 

The liquid argon (LAr) calorimeter surrounds the CTD and is segmented into electromagnetic 
and hadronic sections. It covers the polar angular range 4° < 9 < 154° with full azimuthal 
coverage. The detector is surrounded by an instrumented iron return yoke that is used for muon 
identification (central muon detector CMD, 4° < 9 < 171°). 

The J/ip decay electrons are identified via their energy deposition in the electromagnetic part of 
the calorimeter and via their specific energy loss in the gas of the central drift chambers. Muons 
are identified as minimum ionising particles in the LAr calorimeter or through track segments 
reconstructed in the CMD. 



2.2 Kinematics 

The kinematics for charmonium production are described with the standard variables used for 
deep inelastic interactions, namely the square of the ep centre of mass energy, s = (p + k) 2 , 

'HI uses a right handed coordinate system, the forward {+z) direction, with respect to which the polar angle 9 
is measured, is defined as that of the proton beam. The backward direction (— z) is that of the lepton beam. 
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the squared four momentum transfer Q 2 = —q 2 and the mass of the hadronic final state W = 
\AP + <?) 2 - Here k, p and q are the four-momenta of the incident lepton, proton and virtual 
photon, respectively. In addition, the scaled energy transfer y = p ■ q/p ■ k (energy fraction 
transferred from the lepton to the hadronic final state in the proton rest frame) and the J/ip 
fractional energy z = (p^ ■ p)/(q ■ p) are used, where p^ denotes the J/ ip four-momentum. 

The event kinematics are reconstructed using a method which combines the measurement of 
the scattered lepton and the hadronic final state to obtain good resolution in the entire kinematic 
range. The variable Q 2 = AE E' cos 2 y is reconstructed from the energy E' and angle 9 e of 
the scattered lepton (E is the energy of the incoming lepton). For the calculation of y and z the 
hadronic final state is used in addition. Thus 



Y^jE-p^ p_ ±1 p = {E- Pz ), 

E(E- Pz ) ys/2 Z had (E- Pz y 

where ^2(E — p z ) runs over all the final state particles including the scattered lepton, and in 
Y2had(E ~ on ^y tne fi na l state na drons are summed. The J/ip momentum is reconstructed 
from the momenta of the decay leptons. For the calculation of the sums in equations ([[]) a 
combination of tracks reconstructed in the CTD and energy depositions in the LAr and SpaCal 
calorimeters is used. W is reconstructed using the relation W 2 = ys — Q 2 . 

Differential cross sections are measured as functions of the following variables: Q 2 , W, z, the 
transverse momentum squared of the J/ip with respect to the beam axis p\^ and the rapidity f] 
of the J/ip in the laboratory frame Y\ a b- Differential cross sections are also given for p* 2 ^ and 
Y*, which are computed in the 7*p centre of mass frame. The resolution, as determined from 
the detector simulation, is 2 — 5% for the variables Q 2 , p 2 ^, Y iab and Y*. For z the resolution 
is ~ 8% at high z ~ 1 degrading to 15% at low z values. For W the resolution is ~ 7% for 
W < 145 GeV and ~ 4% above. The resolution of p* 2 ^ is somewhat worse (~ 30% of the 
chosen bin widths). 



2.3 Monte Carlo Simulations 

Corrections for detector effects are applied to the data using a Monte Carlo simulation in which 
the HI detector response is simulated in detail. The simulated events are passed through the 
same reconstruction and analysis chain as the data. The correct description of the data by the 
simulation is checked by independent measurements. Residual differences between data and 
simulation, e.g. in the efficiencies of the lepton identification or of the trigger, are included in 
the systematic uncertainties (Table [I]). 



The Monte Carlo generator used for inelastic J/ip production is EPJPSI [ ]T7| ] which generates 
events according to the Colour Singlet Model in leading order. In contrast to the standard 
version used previously [JTTJ], the full dependence of the matrix element on Q 2 has been imple- 



mented [[T8[]. In order to achieve a good description of the data, the events are reweighted in 



2 The rapidity Y = | In fj^ of the J/ ip is calculated with respect to the proton direction in the laboratory 
frame and with respect to the photon direction in the photon-proton centre of mass frame. 



7 



using a parametrisation of the measured Q 2 distribution. A systematic uncertainty of ±5% is 
estimated for this procedure by repeating the analysis without this reweighting. 

Diffractive production of J/ip and ip(2S) mesons is simulated using DIFFVM [|T9p with pa- 
rameters which have been tuned to HERA measurements. Contributions from the production 
of bb quark pairs with subsequent formation and decay of 6-flavoured hadrons, b — > J/ip + X, 
are simulated by the AROMA Monte Carlo program [pQ|] . The total AROMA cross section is 
normalised to the measured value of 16.2 nb pip. 

2.4 Radiative Corrections 

The measured cross sections are given in the QED Born approximation. The effects of higher 
order processes, mainly initial state radiation, are estimated using the HECTOR program [|2"2|]. 
With the requirement that ^2(E — p z ) > 40 GeV (see below) the radiative corrections amount 
to —(4 ... 5)% and depend only weakly on Q 2 and W. A correction of — (5 ± 4)% is applied. 

3 Data Analysis 
3.1 Event Selection 

Events with Q 2 > 2 GeV 2 are selected by requiring a scattered lepton with a minimum energy 
deposition of 12 GeV in the electromagnetic calorimeter and a lepton scattering angle larger 
than 3°. The z coordinate of the vertex position is determined for each event and required to lie 
in the beam interaction region. In order to minimise the effects of QED radiation in the initial 
state, the difference between the total energy and the total longitudinal momentum {E — p z ) 
reconstructed in the event is required to be larger than 40 GeV. If no particle, in particular no 
radiated photon, has escaped detection in the backward direction, the value of ^2(E — p z ) is 
expected to be twice the incident lepton energy, i.e. 55 GeV. 

The J/ip decay leptons are reconstructed as two oppositely charged particles with transverse 
momenta p t > 0.8 GeV in the CTD. Both tracks have to be identified as muons with polar 
angles in the range 20° < 9 < 160° or as electrons in the range 30° < 9 < 150°. There is a 
considerable non-resonant background, mainly due to misidentified leptons (compare Fig. [|), in 
particular at low values of z. Therefore the number of J/ip candidate events in a given analysis 
interval is extracted by fitting the mass distribution with a superposition of a Gaussian of fixed 
width and position (determined by a fit to all data) to describe the signal and a power law com- 
ponent to describe the background. The number of signal events is then obtained by counting 
the number of lepton pairs in the interval 2.85 < M M/1 < 3.35 GeV and subtracting the fitted 
amount of background in this interval. This method was found to give stable and reliable results 
in most regions of phase space. The statistical error on the number of signal events is estimated 
from the statistical error on the number of events (signal+background) in the mass interval. This 
method leads to a loss of events for the decay of the J/ip to electrons due to radiation of the 
decay electrons in the material of the detector and due to decays J/ip — > e + e _ 7. A correction 
of ~ 10% is applied. A systematic uncertainty of 3 — 7% depending on kinematic variables is 
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Figure 2: Sum of di-muon and di-electron mass spectra for a) Q 2 > 2 GeV 2 and b) Q 2 > 
12 GeV 2 after all selection cuts. The curves are the results of fits of Gaussian distributions for 
the signal and a power law for the non-resonant background. 



estimated for the determination of the signal event numbers by changing the functional form for 
the background. 

After the cuts described above the main background is due to the diffractive production of J / ip 
mesons which is concentrated at high z values. Diffractive J/ip contributions can experimen- 
tally be suppressed in several ways. Previously, inelastic events were selected by requiring the 
hadronic system X, which is produced together with the J/ip meson, to have a high mass [[TT|]. 
In the present analysis, a selection cut z < 0.9 is used to suppress diffractive elastic and proton 
dissociative events. This range corresponds to the region of validity of the theoretical calcula- 
tions which are used for comparison. A further cut is applied, p*^> 1 GeV 2 , where p*^ is the 
transverse momentum of the J/ip in the photon proton centre of mass frame. After this require- 
ment, the background from diffractive J / ip meson production is estimated to be less than 2% 
and is neglected. 



3.2 Contributions from b and i(j(2S) Decays 

After the cuts described in the previous section, the J /if) sample is dominated by 'direct' in- 
elastic J/ip production, in which the J/ip is directly produced from the cc pair in the process 
7* g — ► cc (g). However there remain contributions from both the diffractive and inelastic pro- 
duction of i/j(2S) mesons and the production of b flavoured hadrons with subsequent decays 
to states involving J/ip mesons. Diffractive ip(2S) events are expected mainly at high z val- 
ues while contributions from b — > J/ip + X are expected at low z values. With the cut on 
Pf 2 jp> 1 GeV 2 the remaining background from ip(2S) is mainly due to diffractive events in 
which the proton dissociates and is estimated to be between 6% and 20% in the highest z bin, 
0.75 < z < 0.9, corresponding to 2 — 10% in the total sample. The lower estimate (6%) is 
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based on a Monte Carlo simulation of diffractive ip{2S) production (the simulated contribution 
is shown in Figs. |3]c) and e). Since diffractive ip(2S) production has not been measured in the 
present kinematic range we consider this to be a crude estimate. An analysis of the present data, 
where events with less than five particles are selected as candidates for ip(2S) — ► J/ipn + n~, 
yields an estimate of 20% in the highest z bin. No correction is applied to the data, since the 
dependence on the kinematic variables, in particular on the transverse momentum of the J/ip 
meson, is poorly known. 

The contribution from b — > J/ip + X, which is expected at low values of z, is estimated from 
a Monte Carlo simulation of bb production [ |2~0"| ] using the measured b cross section [ETQ. It is 



estimated to be 17% in the lowest z bin (0.3 < z < 0.45, compare Fig.|]c) corresponding to 
5% in the total sample. Again, no correction is applied to the data due to the poorly known 
dependences on the kinematic parameters. 

Inelastic production of ip(2S) mesons with subsequent decays ip(2S)—> J/ip + X give a fur- 
ther contribution which at present cannot be distinguished experimentally. It is expected to 
contribute over the whole z range and its dependence on the kinematic variables is likely to be 
similar to that of direct J / ip production. It can thus be considered as a normalisation uncertainty. 
In the photoproduction limit this contribution is estimated to be ~ 15% [^p. 

Summarising, the measured cross sections contain in addition to direct inelastic J/ip mesons, 
contributions from diffractive tp(2S) events and b decays which may amount to as much as 17% 
in total. The distributions in the variables studied have not been measured, but are expected to be 
quite different for the two processes and different from those of the direct inelastic J/ip events. 
No correction or systematic error is applied. Inelastic ip(2S) events on the other hand are 
expected to have similar distributions to the inelastic J/ip events themselves; their contribution 
may be of the order of 15% and can be regarded as a normalisation uncertainty. 



4 Results 

Differential cross sections are determined in the kinematic region 2 < Q 2 < 100 GeV 2 
((Q 2 ) = 10.6 GeV 2 ), 50 < W < 225 GeV, p* 2 ^ > 1 GeV 2 and 0.3 < z < 0.9. A second 
set of differential cross sections is determined for a subset with Q 2 > 12 GeV 2 and, in or- 
der to match the Q 2 range, with p 2 ^ > 6.4 GeV 2 . The average value of Q 2 in this sample is 
(Q 2 ) = 30.9 GeV 2 . The distribution of the invariant mass of the two leptons after all selection 
cuts is shown in Fig. ^. The total number of signal events is 458 ± 30 of which 70 ± 1 1 are at 
Q 2 > 12 GeV 2 andp 2 ^ > 6.4 GeV 2 . 

Comparisons between the data and the Monte Carlo simulation (EPJPSI), which is used to 
correct for detector effects, are shown in Fig. ^|. The simulations take into account the two 
lepton proton centre of mass energies according to the luminosity. The simulation is normalised 
to the data in the interval 0.3 < z < 0.9 after reweighting the events in Q 2 and then describes 
all other distributions well. Monte Carlo estimates of contributions from b — > J/ip + X and 
diffractive ip(2S) production are indicated in the z and p* 2 ^ distributions (Fig. |3]c and e). 

The systematic uncertainties in this analysis are typically 15 — 17% and amount to 21% in single 
bins at low z and W. For the double differential cross sections the corresponding error estimate 
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is 21%. The systematic errors are dominated by uncertainties in estimating the number of events 
in regions of high non-resonant background, by the uncertainty in the Monte Carlo calculation 
used for acceptance and efficiency corrections, and by uncertainties in efficiencies for lepton 
identification and triggering. A list is given in Table [I]. 



Source 


Uncertainty 


[%] 


Reconstruction of scattered lepton: angle 




5 


energy 




5 


Reconstruction of decay leptons: track, vertex 




4 


identification 




6 


Number of signal events 


3- 


-15 


Trigger 




7.3 


Monte Carlo model 




5 


Hadronic energy scale 




4 


Radiative corrections 




4 


Integrated luminosity 




1.5 


Decay branching ratio 




2 


Sum 


15- 


-21 



Table 1 : Summary of systematic errors for the single differential J/ tp production cross sections. 
The error on the number of events depends on z and pf^. The total error is the sum of the 
contributions added in quadrature. 



4.1 Differential Cross Sections 

The differential cross sections for inelastic J/ij) production are displayed in Fig. f| as functions 
of Q 2 and p* 2 ^. In Fig. |5] normalised differential cross sections are shown as functions of the 
variables W, z, p 2 ^, p* 2 ^, Yi a b and Y*. The data points are plotted at the mean value of the 
data in each interval. The differential cross sections are also listed in Tables ^| and ^|. The 
results of the calculations by Kniehl and Zwirner [Q], who applied the NRQCD approach to 
electroproduction of J / ip mesons, are shown for comparison. These calculations only include 
2^2 contributions, which is appropriate for z < 0.9. For easier comparison of shapes the 
data and the calculation in Fig. |5| have been normalised to the integrated cross sections in the 
measured range for each distribution. 

The NRQCD calculations shown in the figures include the contributions from the colour octet 
statesf] 3 Si, 3 .Pj=o,i,2, X Sq as well as from the colour singlet state 3 Si (labelled "CO+CS"). The 
contribution of the colour singlet state is also shown separately (labelled "CS"). The calculations 
depend on a number of parameters. The values used for the non-perturbative long range transi- 
tion matrix elements (LDMEs) were determined from the distribution of transverse momenta of 

Spectroscopic notation is used: 2S+1 Lj where S, L and J denote the spin, orbital and total angular momenta 
of the cc system. 
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J / ip mesons produced in pp collisions [[7J]p]. The bands in Figs. |], |] and |^ indicate the uncertainty 
in the theoretical calculation [^5[]. They cover a charm quark mass of m c = 1.5 ±0.1 GeV, vari- 
ation of renormalisation and factorisation scales by factors 0.5 and 2, the errors of the LDMEs 
as well as the case that either of the two parts of Mr (see footnote H) does not contribute. 



Furthermore, the effect of using the CTEQ5M [24] set of parton density functions instead of 



MRST98LO [23] is included. 



The colour octet contribution dominates the predicted cross section for all values of Q 2 and p* 2 ^ 
(Fig. |]a and c). In order to facilitate the comparison with the data, the ratio data/theory is shown 
on a linear scale in Fig. |]b and d together with a band indicating the uncertainty in the NRQCD 
calculation with CO+CS contributions. The NRQCD calculation overshoots the data by about 
a factor of 2 at low Q 2 , which is at the limit of the theoretical and experimental error. The 
agreement between the data and the theory improves towards higher Q 2 where the theoretical 
uncertainties diminish. For p* 2 ^ similar agreement between data and NRQCD calculation is 
observed. Compared with the colour singlet contribution alone, the data exceed the calculations 
by a factor ~ 2.7, approximately independent of Q 2 , while for pf^ the ratio increases towards 
higher values of pfi. 

In Fig. |5] the measured and the theoretical differential cross sections are normalised to the in- 
tegrated cross sections in the measured range for each distribution. The W and the distri- 
butions (Fig. |]b and f) are reasonably well described in shape by the full NRQCD calculation 
and also by the colour singlet contribution alone, whereas neither fully describes the Y* dis- 
tribution. The z distribution is very poorly described by the full calculation including colour 
octet contributions, while the colour singlet contribution alone reproduces the shape of the data 
rather well. A similar discrepancy between data and NRQCD calculations was observed at high 
z values in the photoproduction limit [|9|,[T(],|2(J]. It may be due to phase space limitations at 



high z for the emission of soft gluons in the transition from the colour octet cc pair to the J/ip 
meson, which are not taken into account in the calculation. In photoproduction, the rapid rise of 
the colour octet contributions towards high z values was shown to be damped after resummation 
of the NRQCD expansion [|2"6|, j2"7[] . The shapes of the p 2 ^ and pf^ distributions (Fig. |5]c and 



e) are rather well described by including CO+CS contributions while the CS contribution alone 
decreases too rapidly towards high values of p 2 ^ or p* 2 ^. Note, however, that higher orders are 
expected to contribute significantly at high values of p t ^ as observed in next-to-leading order 
CSM calculations in the photoproduction limit [@]. 

At higher Q 2 values the theoretical uncertainties of the NRQCD calculation decrease (see 
Fig. |]b). It is therefore interesting to compare data and theory at higher Q 2 . The results for 
Q 2 > 12 GeV 2 (with p 2 ^ > 6.4 GeV 2 ) are given in Fig. § and Table |[ The requirement 
pfy> 1 GeV 2 is retained. The average (Q 2 ) = 30.9 GeV 2 is larger than the squared mass 
of the J /if) meson, ni 2 ^. The statistical precision of these data is limited and no substantial 
change in the comparison of data and theory is seen compared to Fig. £j 

4 The extracted values for the LDMEs depend on the parton density distributions. For the set MRST98LO [ |23| ] 
the values are, in the notation of [|], {O j ^'[ 3 s[ 1} }) = 1.3 ±0.1 GeV 3 , (O j W[ 3 S[ &) ]) = (4.4 ± 0.7) • 10~ 3 GeV 3 
and M^' = (8.7 ± 0.9) • 10~ 2 GeV 3 . is a linear combination of two elements, M r J/V ' = (O j ^ [ x 5^ 8) ]) + 

r(O J /*[ 3 P (8) ])/m2, with r = 3.4. 
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4.2 Double Differential Cross Sections 



In the calculations the relative contributions of the colour octet states to the cross sections vary 
with z, Q 2 and pft (compare Figs. 0, |5|a and c). Therefore, differential cross sections do jdp* 2 ^ 
and da/dQ 2 are determined^ in three intervals of z and compared with that for the whole z 
range in Fig. [7| The dependence on Q 2 and p* 2 ^ is seen to be similar in the three z regions. In 
order to make a quantitive comparison, the differential cross sections for the whole z range are 
fitted with functions oc (Q 2 + m 2 ,^~ n or oc (p*^ + m 2 ^)'™ yielding (n = 3.36 ± 0.53) and 
(m = 4.15 ± 0.50), respectively, where total experimental errors are given. The results of these 
same fits are then compared with the data in the three z intervals after normalising the curves 
at low Q 2 or p* 2 ^. In Figs. [7]b and d the ratio of the data over the scaled fit is shown. The data 
in the three z bins are reasonably described by the same functional form although there is an 
indication of a faster fall with Q 2 at high z than in the total z range. In view of the contributions 
at high z from diffractive ip(2S) production, which are expected to have a different dependence 
on Q 2 , firm conclusions cannot be drawn. The observed dependence on p* 2 , is within errors the 
same as that observed in photoproduction (m ps 4.6 ± 0.1) p5[|. 



4.3 7*p Cross Sections and Comparison to Photoproduction 

For comparison with results in the photoproduction limit the cross section for 7*p — > J/ip X as 
a function of W is calculated by dividing the ep cross section by the photon flux integrated over 
the analysis intervals [|2"g[]. 

The total cross section cr{^*p — > J/ipX) is shown as a function of W in Fig. [8| and listed 
in Table ^ for the present data ((Q 2 ) = 10.6 GeV 2 ). It is compared with the cross section in 
the photoproduction limit ((Q 2 ) ~ 0.05 GeV 2 ) in an otherwise similar kinematic range [|26|]. 
Parametrising the cross section as (W/lVo)' 5 yields 5 = 0.65 ± 0.25 for the present data, where 
the total experimental error was used in the fit. The value is consistent with that obtained in 
photoproduction (0.49 ± 0.16 [E^p). The W dependences are expected to be similar since they 



reflect the x dependence of the gluon distribution with a scale ~ Q 2 + m 2 



4.4 Decay Angular Distributions 

Measuring the polarisation of the produced J / ip meson has been proposed as a means of dis- 
tinguishing the various CO and CS contributions to J/tp production [^],|Sp. The polar (9*) decay 
angular distributions are measured in the rest frame of the J/tp using the J/ip direction in the 
7*p system as reference axis (helicity frame). They are shown in Fig. |9| (and listed in Table [7]) 
for the whole Q 2 range and separately for regions of low and high Q 2 . The cos^* distribution 
is expected to have the form 

< ifT oc 1 + a cos 2 9*. (2) 



d cos 9* 



5 The corresponding double differential cross sections da / dpfydz and da/dQ 2 dz are listed in Table |5|. 
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A value of |a| £ 0.5 is expected, where a can be negative, zero or positive depending on 
which intermediate cc state dominates the production [Q] . Fitting the data distributions with a 
function of the form (Q) yields a value of a = — 0.1+q 3 m the whole Q 2 range (Fig. ^|a). For 
2 < Q 2 < 6.5 GeV 2 (Fig. gb), a = -0.4^j is found and for 6.5 < Q 2 < 100 GeV 2 (Fig. |c) 
a = 0.7^Q g. The total experimental errors were used in the fits. Although the central values 
suggest a change from a negative to a positive value of a as Q 2 increases, this tendency is not 
significant. Predictions using the k t factorisation approach flTSp, shown in Fig. are compatible 
with the measurements. 



5 Summary and Conclusions 

A new analysis of inelastic electroproduction of J /if) mesons has been presented. Due to the 
increased statistics the kinematic range has been extended to 50 < W < 225 GeV and reaches 
average values of Q 2 larger than the squared mass of the J/ip meson. The cross sections are 
measured in the range 0.3 < z < 0.9 and p^> 1 GeV 2 where direct inelastic J/ip production 
dominates. Differential cross sections at average values (Q 2 ) = 10.6 and 30.9 GeV 2 are pre- 
sented as functions of Q 2 , W, z,p 2 ^, pf^, Y and Y*. Recent theoretical calculations by Kniehl 
and Zwirner [Q] in the framework of the non-relativistic QCD (NRQCD) approach including 
colour octet and colour singlet contributions ("2 — ► 2" diagrams) are compared with the data. 
At both average Q 2 values, reasonable agreement is observed in the shapes of most distributions 
except that of z, which is described much better by the colour singlet contribution alone (in a 
recent resummation of soft gluon processes a similar observation in the photoproduction limit 
could be explained through a damping of the rapid rise of the colour octet contributions towards 
high z values). The absolute value of the full NRQCD cross section is a factor ~ 2 above the 
data at low Q 2 but approaches the data at higher Q 2 to within 15% which is well within exper- 
imental and theoretical uncertainties. The colour singlet contribution alone is roughly a factor 
2.7 lower than the data. The differential cross sections in p 2 ^ and p* 2 ^ are described better 
when CO contributions are included. In the photoproduction limit a successful description of 
the p\^ spectrum has been achieved within the Colour Singlet Model by including NLO cor- 
rections. These corrections are, however, not yet available for the electroproduction case under 
consideration here. The dependence of the YP cross section on W is the same, within errors, 
as in the photoproduction case. 

Further distributions are studied in an attempt to assess the relative importance of the different 
CO and CS terms. Since their contributions are expected to vary with z, differential cross 
sections da/dQ 2 and da/dpf^ are measured in intervals of z. The shapes of the p* 2 ^ and Q 2 
spectra are found to be similar to those over the whole z range although there is an indication of 
a steeper Q 2 dependence at high z. A fit to the distribution of the polar decay angle in the helicity 
frame covering the whole Q 2 range yields a = — O.ligi f° r a parametrisation 1 + a cos 2 9*. 

Altogether the measurements presented here provide significant new information which will 
aid the further development of a quantitative understanding of J/ ip meson production within 
pQCD. 
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Figure 3: Comparison between data and Monte Carlo simulations for ep — > e J ftp X in the 
region 50 < W < 225 GeV, p*J > 1 GeV 2 and 0.3 < z < 0.9 after all selection cuts and 
after subtraction of non-resonant background. Distributions of a) Q 2 , b) W, c) z, d) p 2 ^, e) 
pl 2 ^, f) F* and g) F; a 5 are shown. The error bars on the data points are statistical only. The full 
histograms show the inelastic J/ ip events from EPJPSI after weighting the events according to 
a parametrisation of the Q 2 dependence of the data and normalising over the z range measured. 
The estimated contribution from b decays (AROMA, dashed) is also shown in c) and from 
diffractive V ; (2S')(EPJPSI, dotted) in c) and e). These contributions are not included in the full 
histograms. 
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Figure 4: Differential cross sections a) da/dQ 2 and c) da /dp* 2 ^ for the inelastic process ep 



e J X in the region 50 < W < 225 GeV, Q 2 > 2 GeV 2 , pjj, > 1 GeV 2 and 0.3 < 2 < 0.9. 
The NRQCD calculation is shown for comparison (CO+CS, light band) and the colour singlet 
contribution separately (CS, dark band). In b) and d) the ratio of data/theory is shown for the 
two cases. The theoretical uncertainty in the full calculation (CS+CO) is shown as a band 
around 1. The inner error bars of the data are statistical, the outer error bars contain statistical 
and systematic uncertainties added in quadrature. 
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Figure 5: Normalised differential cross sections for the inelastic process ep — ► e J/^ X in 
the kinematic region 2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, p*J > 1 GeV 2 and 0.3 < 
z < 0.9. a) l/ada/dz, b) 1/ada/dW, c) l/ada/dp^, d)l/a da/dY* e) 1/ada/dp^ and 
f) 1 / a da / dY iab . The inner error bars are statistical, the outer error bars contain statistical and 
systematic uncertainties added in quadrature. The histograms show calculations for inelastic 
J/ip production within the NRQCD factorisation approach ^\ which have been normalised to 
the integrated cross section. The light band represents the sum of CS and CO contributions 
and the dark band the CS contribution alone (both are separately normalised). The error bands 
reflect the theoretical uncertainties (see text). 
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Figure 6: Normalised differential cross sections for the inelastic process ep — > e J/^ X in 
the kinematic region 12 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, p 2 ^ > 6.4 GeV 2 , p*J > 
1 GeV 2 and 0.3 < z < 0.9. a) l/ada/dz, b) 1/ct Ax/cW, c) 1 / 'a da / 'dp 2 t ^, d) 1/adajdY* 
e) l/ada/dp 7 ^ and f) 1/ada/dQ 2 . The inner error bars of the data points are statistical, 
the outer error bars contain statistical and systematic uncertainties added in quadrature. The 
histograms show calculations for inelastic J/ij) production within the NRQCD factorisation 
approach [§]. The light band represents the sum of CS and CO contributions and the dark 
band the CS contribution alone (both are separately normalised). The error bands reflect the 
theoretical uncertainties (see text). 
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Figure 7: Differential cross sections for ep — > e J/ipX in three z intervals and in the full z 
range, a) da/dpf^ and c) da/dQ 2 for low (0.3 < z < 0.6, open points), medium (0.6 < z < 
0.75, triangles) and high (0.75 < z < 0.9, squares) values of z in comparison with the results 
for the full z region (full points). The inner error bars indicate the statistical uncertainty, while 
the outer error bars show the statistical and systematic uncertainties added in quadrature. For 
clarity, the data have been scaled by the factors indicated. The data in the complete z range are 



parametrised by fits of the form (Q 



-m 



J/i>> 



~ n and (p*f^+m 



\ The same parametrisations 



are also shown for the data in the three z bins after normalising them to the data at low Q 2 or 
■p* 2 ^. In b) and d) the ratio of the data to these parametrisations is shown on a linear scale using 
the same symbols as in a) and c). Note that for clarity the points have been shifted in Q 2 and 

n* 2 
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Figure 8: Total cross section for 7* p -> J/^Xat(Q 2 ) = 10.6 GeV 2 in the range 0.3 < z < 0.9 
as a function of W. Photoproduction data with similar cuts in z and p 2 ^ are included for 
comparison. The curves are the results of fits with the function (W/W ) s . The inner error bars 
on the points indicate the statistical uncertainty, while the outer error bars show the statistical 
and systematic uncertainties added in quadrature. 
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Figure 9: Differential cross sections 1/ a da / dcosO* in ep — > e J/ip X in the kinematic region 
50 < W < 225 GeV, p*J > 1 GeV 2 and 0.3 < z < 0.9 normalised for | cos0*| < 0.9. 
a) 2 < Q 2 < 100 GeV 2 , b) 2 < Q 2 < 6.5 GeV 2 , c) 6.5 < Q 2 < 100 GeV 2 . The inner error 
bars indicate the statistical uncertainty, while the outer error bars include the statistical and 
systematic uncertainties added in quadrature. The shaded regions show the result of fits with 
the form ~ 1 + a cos 2 9* and correspond to a variation of the fit parameter a by ±1 standard 
deviation. The dashed lines are the result of a prediction using the k t factorisation approach jJTj]. 
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Table 2: Differential cross sections with statistical and systematic errors in the range 
2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 0.3 < z < 0.9 and p* t % > 1 GeV 2 . 
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Table 3: Differential cross sections with statistical and systematic errors in the range 
2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 0.3 < z < 0.9 and p* t % > 1 GeV 2 . 
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Table 4: Differential cross sections with statistical and systematic errors in the range 
12 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, p 2 ti , > 6.4GeV 2 , 0.3 < z < 0.9 and 
j^>lGeV 2 . 
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Table 5: Double differential cross sections with statistical and systematic errors in the range 
2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 0.3 < z < 0.9 andp* 2 ^ > 1 GeV 2 . 
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Table 6: Cross sections with statistical and systematic errors for 7* ± p — > J/i/) + X at 
(Q 2 ) ~ 10.6 GeV 2 in the range 0.3 < z < 0.9 and p* t % > 1 GeV 2 . 



e + p -»■ e + J/<0 + X 
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125. ±20. ±20. 


0.55-0.90 


0.70 
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0.20-0.45 


0.32 
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6.5 < Q 2 < 100 GeV 2 
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-0.66 


52.4 ± 11.0 ±8.4 
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Table 7: Differential cross sections with statistical and systematic errors in the range 
2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 0.3 < z < 0.9 andp* 2 ^ > 1 GeV 2 . 
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